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ABSTRACT 

Dwarf novae in globular clusters seem to be rare with only 13 detections in the 157 known Galactic 
globular clusters. We report the identification of a new dwarf nova in M13, the 14th dwarf nova 
identified in a globular cluster to date. Using the 2m Faulkes Telescope North, we conducted a search 
for stars in M13 that show variability over a year (2005-2006) on timescales of days and months. 
This led to the detection of one dwarf nova showing several outbursts. A Chandra X-ray source 
is coincident with this dwarf nova and shows both a spectrum and variability consistent with that 
expected from a dwarf nova, thus supporting the identification. We searched for a counterpart in 
Hubble Space Telescope ACS/WFC archived images and found at least 11 candidates, of which we 
could characterize only the 7 brightest, including one with a 3a Ra excess and a faint blue star. The 
detection of one dwarf nova when more could have been expected likely indicates that our knowledge 
of the global Galactic population of cataclysmic variables is too limited. The proportion of dwarf 
novae may be lower than found in catalogs, or they may have a much smaller mean duty cycle (~1%) 
as proposed by some population synthesis models and recent observations in the field. 
Subject headings: Galaxy: globular clusters: individual: M13 (NGC 6205) - X-rays: general - Stars: 
dwarf novae - Stars: novae, cataclysmic variables 



1. INTRODUCTION 

Globular clusters (GCs) are old, gravitationally bound 
stellar systems with extremely high stellar densities in 
their core which enable cluster members to dynamically 
interact through regular encounters. Close binary stars 
have been suggested to pl ay a significant r ole in the dy- 
namical evolution of GCs (|Hut et al.| [l992). Encounters 
between passing stars and close binaries may harden the 
latter and accelerate the passing star. This leads to the 
formation of a variety of tight and interacting binaries, 
and by accelerating stars in the core, close binaries may 
delay or halt the inevitable gravitational core collapse of 
GCs. 

Many of the X-ray sources in GCs discovered in the 
70's have been resolved into a multitude of close bina- 
ries and relat ed systems with XM M-Newton and Chan- 
dra (see e.g. iServillat et a l. 2008a.b). These include ac- 
tive and quiescent low-mass X-ray binaries (LMXBs) 
with a neutron star (NS) primary, millisecond pulsars 
(MSPs), cataclysmic variables (CVs) and magnetically 
active binaries. The number of NS LMXBs in GCs 
as well as part of the CV population appears to scale 
with the stellar encounter rate of their host cluster, in- 
dicating that the y are indeed produced through dynam- 
ical mechanisms (jGendre et al.l 120031 : iPooley et al .1120031 : 
iPoolev fc Hutl 120061 ). The same correlation is observed 
with MSPs in GCs detected in radio (jHui et al.lfcOlOf ). 
and even in 7-rays (jAbdo et al.ll2010f ). 

CVs are semi-detached binary stars with orbital peri- 



ods typically of the order of hours, consisting of a white 
dwarf primary accreting via Roche lobe overflow from a 
companion which is usually a late- type, generally main- 
sequence star. Subtypes of CVs include non-magnetic 
CVs where an accretion disc forms around the primary, 
and magnetic systems where the white dwarf magnetic 
field is strong enough to disrupt (polar CVs, B > 10 T G) 
or truncate (intermediate polars, 10 5 < B < 10 7 G) the 
accretion disc, and the mass fl ow is then channeled along 
the m agnetic field lines (e.g. iCropperl 119901 : iPattersonl 
1994). For systems with an a ccretion disc , it is believed 
that a thermal instability (e.g. Osaki 1996h is the cause of 
repetitive outbursts observed in CVs called dwarf novae 
(DNe). Optically, DN eruptions have amplitudes of 2-6 
magnitudes in V, a duration of a few to 2 days and a 
recurrence timescale of weeks to years (e.g. I0sakil ll996). 

From the Catalog and Atlas of Cataclysmic Variables 
(jDownes et al J I2QQ1L updated 2006), DNe are the most 
common CVs identified in the field, with about half of 
the identified CVs being DNe. However, the sample of 
known field CVs suff ers from various selection effects 
(jPretorius et al. 2007), and is biased towards bright and 
highly variable systems, such as DNe. We note that 
the distribution of those CVs is almost isotropic, which 
would indicate that the known population is mainly local 
(less than 1 kpc around the Sun). The global Galactic 
population and its intrinsic composition is thus not fully 
understood. 

Population synthesis models indicate that a fainter 
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population showing fewer outb ursts may dominate the 
popu l ation of DNe (up to 70%, lKolblll993b iHowell et all 
1997; iWillems et all l2QQ5f ). Recent detections in the 
SDSS of faint CVs with short periods showed one of 
the first hints for th e existence of such a population 
(|Gansicke et al.l l2009). but its proportion is still not con- 
strained. 

If we limit our knowledge to the locally known popula- 
tion of CVs, t here appear to be f e wer DN outbursts from 
CVs in GCs (jShara et all 11993 iDobrotka etHI 12006). 
This would mean that the population of CVs in GCs 
is intrinsically different from the known population in 
the field. However, it is clear that the known popu- 
latio n in the field is no t a reliable reference. We note 
that Sha ra et all (| 19961 ) based their comparison on a 
small sample of about 20 well studied DNe which are 
not likely to represent the global population. Two hy- 
potheses remain to be fully tested, (i) Dynamical ef- 
fects and a specific environment in GCs (low metallic- 
ity, long timescales) in GCs could shape differently the 
population of CVs. For example, it has been proposed 
that CVs in GCs are dominated by moderately mag- 
netic systems (intermediate polar s) with low accretion 
rates which sho w fewer outbursts (Do brotka et al.l 12006; 
lGrindlavl[T999h . Also, typical masses of the accreting 
primaries may be higher, as massive white dwarfs can be 
more efficiently pl aced in binaries by tidal capture and 
3-body processes (|Ivanova et al.l 12006). (ii) The lack of 
DN outbursts from CVs in GCs could be due to detection 
biases. Due to the crowding in optica l, DN outbursts are 
indee d more difficult to catch in GCs (jPietrukowicz et all 
2008). These questions can only be solved with further 
observations and comparisons with more complete pop- 
ulations of CVs, both in the field and in GCs. 

We present here a search for optically variable stars, 
and particularly DN outburs ts, in t he GC M13 (NGC 
6205). From the catalog of Harris (1996, updated in 
2010), M13 is at a distance of 7.1 kpc from the Sun. 
Its core radius is 0. / 62, which corresponds to the half 
width half-maximum of the luminosity profile of the clus- 
ter, and its half-mass radius is 1.'69 (region encircling 
half of the mass, approximated to be half of the light 
from the cluster). M13 is a low metallicity cluster with 
[Fe/H] = —1.53. The reddening along the line of sight is 
E(B — V) = 0.02, which corresponds to a low extinction 
Ay ~ 0.06 and an absorption from a column density of 
N H - 1.1 x 10 20 cm 2 . 

We first present our search for optical variables and 
outbursts (©. We then detail Chandra X-ray observa- 
tions of M13 used to search for X-ray sources coincident 
with the optically variable ones (§2). In $4j we describe 
one peculiar optical event aligned with an X-ray source 
and present archived Hubble Space Telescope data to un- 
cover possible counterparts. Finally, we discuss the re- 
sults in $5l 

2. OPTICAL OBSERVATIONS 

2.1. Data and reduction 

The GC M13 has been observed with the Faulkes Tele- 
scope North (FTN) at Haleakala on Maui, USA, dur- 
ing the nights indicated in Table [TJ The FTN is a 2m 
robotic telescope optimized for research and education 
(jRoberts et a l. 2008; Lewis fc R oche 20091) and is part of 



Table 1 

Faulkes Telescope North U-band 
observations. 



Epoch 


Date 


MJD 


Exposures 


1 


2005-02-18 


53419 


60s, 120s 




2005-02-22 


53423 


120s 


2 


2005-03-02 


53431 


3x120s 




2005-03-04 


53433 


120s 


3 


2005-04-12 


53472 


8x120s 




2005-04-15 


53475 


9x120s 




2005-04-16 


53476 


9x120s 




2005-04-17 


53477 


9x120s 


4 


2005-05-05 


53495 


2x800s 


5 


2006-03-12 


53806 


44x200s 



the Las Cumbres Observatory Global Telescope Network. 

Imaging was obtained in U, B, V, R, i, Ha and Oil I 
filters between February 2005 and April 2006. We fo- 
cus here on the U band as this dataset is the cleanest 
and the most complete, and CVs are generally blue ob- 
jects so they are likely to be among the brightest objects 
compared to GC main sequence stars in the U band. 
Moreover, narrow band filters are not as deep, and we 
obtained fewer images in the other filters. 

The DillCam CCD detector we used has 2048 x 2048 
pixels binned 2x2 prior to readout into effectively 1024 x 
1024 pixels. The field of view is 4.7 x 4.7 and the pixel 
scale is 0'/278 per pixel. Automatic pipelines de-bias and 
Hatfield the Faulkes science images using calibration files 
from the beginning and the end of each night. 

The seeing ranged from V!2 to more than 3". We 
discarded images with poor signal-to- noise ratio (S/N; 
mostly due to thin clouds) and images with tracking 
problems (where the stars are non-spherical and stretch 
more than a few pixels during the exposure), an excess 
of hot pixels, or a background gradient across the im- 
age. Finally, 78 U-band images were kept. They can be 
grouped into 5 epochs of observations, each epoch last- 
ing one to a few days. The first 4 epochs were taken over 
two and a half months and the fifth epoch about one year 
later. 

2.2. Search for variability 

We looked for variable stars in Ml 3 using difference 
imaging photometry. First, we aligned the images by 
selecting 3 stars in a reference image. The task imregister 
under IRAFQ was used to get a set of 1024 x 1024 pixels 
images all aligned on the reference image after rotation 
and translation. Due to this method, the edges of some 
images were cut, and the usable field of view for this 
search has been reduced to ~ 4' x 3.'7, which still included 
the half-mass radius of M13. 

We computed a reference image by taking the median 
of 5 images with the best seeing (about l'/2). This im- 
age was then subtra cted from all the other images using 
the ISIS 2 software (|Alardll200tih. This tool uses adap- 
tive kernel techniques ( Alard fc Luptonlll998f ) and gives 
optimal results in the case of crowded fields, variable see- 
ing and point spread function (PSF). We set the kernel 
size to 9 x 9 and the stamp size to 15 x 15 pixels (the 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation 
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Figure 1. Lightcurves of 3 variable star candidates associated 
with BL Her variables, see Table [2] The Y axis gives the flux 
relative to the reference image in arbitrary units with one sigma 
error bars. 
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Figure 2. Lightcurves of 3 variable star candidates associated 
with RR Lyr stars, same comments as in Figure ^ 

stamp size is the area taken by the program around each 
object for computation). Three Gaussians were used in 
the convolution. The background is nearly flat so it was 
fitted with a constant. The spatial variations in the PSF 
kernel were modeled by a 4th order polynomial. After 
this process, the mean of the obtained images was close 
to zero, as expected, with some isolated positive or neg- 
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Figure 3. Lightcurves of 3 variable star candidates associated 
with RG stars, same comments as in Figure ^ 
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Figure 4. Lightcurves of 3 variable star candidates not associated 
with previously known variable stars, see Figure^ Previously, star 
7 had only been suspected to be variable. 

ative sources. 

For each subtracted image, we detected variable star 
candidates using IDL procedures from the AstrolitB (e.g. 
find). We estimated the background level and the stan- 
dard deviation a in an annulus around the source and 

2 http://idlastro.gsfc.nasa.gov/ 
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Figu re 5. M 1 3 obs erved with the FTN (reference image). The core and half-mass radii are shown in black. Known variable objects from 
Kopacki et al. (2003) are plotted as blue squares, X-ray sources detected with Chandra as red triangles and variable stars detected in this 
work as black circles. 



performed a detection for sources with a significance 
higher than 5<j. However, when irregular features were 
present in the subtracted images, the IDL/Astrolib find 
procedure gave a high number of false detections. When 
this was the case, all the detections for that image were 
ignored. We thus selected the most secure detections 
and possibly missed low significance events. We finally 
merged the different lists and found a total of 16 variable 
star candidates in the field of view. They are listed in 
Table [2] 

For each detected variable star in each subtracted im- 



age, we first fitted the star with a Gaussian, and esti- 
mated the integrated flux in arbitrary units. The typical 
observed full width half-maximum (FWHM) value in the 
subtracted images is about 1//5, slightly larger than the 
FWHM of the reference image. We estimated the signif- 
icance of the variable star candidates by comparing their 
amplitude with the standard deviation of the surround- 
ing background. For the light curve extraction, we kept 
only measurements that were higher than three times 
the standard deviation for each variable star. Some ex- 
amples of lightcurves are shown in Figures [TJ [2j [3l and 
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Table 2 

List of variable stars. 



# 


1 
i 


RA 
[h:m:s 


1 
J 


L 


Dec 
d:m:s] 




Name 


Type 


Period 
fdavl 


<V> 


AV 


1 


16: 


■A 1 


,ou. 


.90 


36: 


:27:48. 


.7 


v02 


BL Her 


5.11 


13.01 


0.87 


2 


16: 


:41 


:44. 


.57 


36: 


28:04. 


4 












3 


16: 


:41 


;42. 


.73 


36: 


27:31. 


1 


v25 


RR Lyr 


0.43 


- 


0.44 


4 


16: 


:41 


;42. 


.46 


36: 


:28:08. 


1 




X-ray 








5 


16: 


:41 


:41. 


.50 


36: 


:27:47. 


,8 


v35 


RR Lyr 


0.32 




0.25 


6 


16: 


:41 


:36. 


.62 


36: 


26:35. 


2 


vll 


RG 


91.77 


11.93 


0.13 


7 


16: 


:41 


:51. 


.99 


36: 


26:29. 


,9 


c 


U-V 








8 


16: 


:41 


:38. 


.68 


36: 


25:37. 


.2 


v38 


RG 




12.12 


0.07 


9 


16: 


:41 


:48. 


.00 


36: 


29:10. 


1 


v06 


BL Her 


2.11 


14.08 


0.60 


10 


16: 


:41 


:50. 


,95 


36: 


:28:54. 


,6 


vl7 


RG 


43.04 


11.98 


0.38 


11 


16: 


:41 


:46. 


,36 


36: 


:27:39. 


.9 


v05 


RR Lyr 


0.38 


14.77 


0.42 


12 


16: 


:41 


:46. 


,49 


36: 


:27:28. 


.1 


vOl 


BL Her 


1.46 


14.09 


1.04 


13 


16: 


:41 


:49, 


,07 


36: 


:27:08. 


.0 


v34 


RR Lyr 


0.39 


14.83 


0.32 


14 


16: 


:41 


:37 


,11 


36: 


:26:28. 


.7 


v07 


RR Lyr 


0.31 


14.93 


0.30 


15 


16: 


:41 


:35. 


,51 


36: 


:27:27. 


,6 


v42 


RG 




11.94 


0.10 


16 


16: 


:41 


:46. 


,29 


36: 


:27:37. 


.9 


v09 


RR Lyr 


0.32 


14.82 


0.43 



Note. — Variable stars detected with their coordinates (columns 2 and 
3). When availa ble, we list the nam e in column 4, and properties of the 
stars derived by [Kopacki et al. ( 2003) in columns 5 to 8: the type of star, 
the period, the average Vega magnitude in the V band, and the amplitude 
of the variabi l ity A V observed. St ar 7 is labelled "c" as first reported 
bv IMein unger (1978) and reused by Kopacki et al. (2003 ). RG stands fo r 
red giant and U-V for U-V-bright object as studied by|Zinn et al. (|19720 . 
X-ray indicates an association with an X-ray source. RR Lyr are short- 
period (0.2 to 2 days) population II pulsating variables often found in GCs. 
BL Her are a subclass of Type II Cepheid variables with periods between 
1 and 8 days. 



[H The zero level in those plots corresponds to the flux 
of the star in the reference image. If no point is reported 
for a given day, it is either because the luminosity of 
the object is at the same level as in the reference image 
(then no emission is seen in the subtracted image), or the 
emission is not significant (< 3<r) compared to the sur- 
rounding background of the subtracted image. We thus 
have discarded some low significance measurements. 

In order to check if our method could lead to false 
positives, we added manually 5 entries in the variable 
star candidate list (aligned with bright stars or randomly 
placed in the image). The chosen positions did not cor- 
respond to objects detected as variables. We applied the 
same extraction process for all 5 positions and we did 
not obtain any significant measurements, i.e. we did not 
obtain false detections. 

2.3. Known variable stars 

In order to compare our list of variable stars with previ- 
ous catalogs, we added an astrometric solution to the ref- 
erence image by aligni ng stars using the 2MASS Catalog 
(|Skrutskie et al. 2006) which provided 15 good reference 
stars. We used GAIA0 to fit the position of those stars 
and generate an astrometric solution for the reference 
image with an RMS error of 0'/2. Taking into account 
the accuracy of 2MASS, the total astrometry over the 
field of view is 0'/3 (1<j error). 

We compared our variable star list t o the Catalog 
of Variable Stars in Glob ular Clusters ([Clement et al I 
120011 : iKopacki et all 120031 ). We note that, as the sub- 

3 GAIA was created by the now closed Starlink UK project and 
has been more recently supported by the Joint Astronomy Cen- 
tre Hawaii funded again by the Particle Physics and Astronomy 
Research Council (PPARC) and more recently by its successor or- 
ganisation the Science and Technology Facilities Council (STFC). 



tracted images are convolved (and thus smoothed) during 
the subtraction process, the positional accuracy can be 
lower than for the reference image. We iteratively cross- 
matched the two lists to find a systematic offset. This 
bore-sight correction amounts to 0'/93 and 0'/66 in RA 
and Dec, respectively, with a la error of 0"15 between 
the 2 sets, which adds up to the astrometric error. This 
led to the identification of 14 variable star candidates as 
already known variables. The name, type, period and 
magnitude of the known variable stars are reported in 
Table [21 The 16 variable stars detected in this study as 
well as the known variables stars and the X-ray sources 
in the field of view are shown in Figure [5] 

Although our data are insufficient to determine a pe- 
riod, the variability observed for stars 1 (v02) and 9 (v06) 
presented in Figu re [1] is consistent w ith the previously 
observed periods (Ko packi et al. 2003). The same trend 
is seen for star 3 (v25) in Figure [1] with a shorter period 
of 0.3 day and all the other lightcurves extracted. 

Red giants (RG) stars have pre viously been observe d 
with semi- or irregular variability (Kopacki et al. 2003). 
The lightcurves we obtained for this category of candi- 
dates are consistent with this observation (Figure [3]) but 
due to our limited sampling they could as well be peri- 
odic. 

The variable star 7 is a known U and V bright ob - 
ject (misleadingly called UV-bright bv lZinn et al.lll972[ ), 
associated with M13 based on radial velocit y measure- 
ments. The U and V bright stars referred bv IZinn et al.l 
(|1972[ ) are in fact related to post-horizontal branch stars 
when plotted in a CMD (see their Figure 2). Star 
7 corresponds to L99 3 and IV- 5 2 in the catalogs of 
iLudendor fi (1905) and I Arpl JT955), respectively and has 
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Figure 6. Chandra image of M13 with circles corresponding to 
the core (0. / 62) and half-mass radii (L'69). Sources are labelled 
according to their distance to the center of the cluster. Labelled 
circles correspond to the 3cr errors including the absolute pointing 
accuracy. 



Table 3 

Detected X-ray sources and expected 
background sources in M13. 



Radius 


0. / 62 


1.'69 


4.' 2 


Detected 


8 


15 


35 


Expected 


0.3 ±0.3 


2.5 ±1.6 


15.4 ±3.9 


Difference 


7.3 ±0.3 


12.5 ±1.6 


19.6 ±3.9 



a spectral type A5III-IV from SIMBA lfl This star has 
been reported as possibly variable (Mei nungerl 1978) 
bu t needed confirmation . No variability was reported 
by iKopacki et al.l (|2QQ3l star labelled c), but we show 
its variability in this study, which seems to occur on 
timescales of a few days. 

Finally, the stars 2 and 4 have no previous detections 
as variable stars in the optical. Those stars are unlikely 
to be periodic variable stars with periods 0.2 to 8 days as 
previous, more complete searches for such pe riodic stars 
did not detect them (e.g. IKopacki et aLll2003[ ). However, 
we note that the star 2 is variable over short timescales 
of one or two days, which seems quite short compared 
to typical DN outbursts. In order to constrain the na- 
ture of those new events, we searched for possible X-ray 
counterparts, as reported in the next section. 

3. Chandra X-RAY OBSERVATION OF M13 

In the Chandra X-ray Observatory archive, we found 
two observations of M13 with the ACIS-S detector on 
2006 March 10 and 11 (MJD 53804-53805) for 27 and 
28 ks with a gap of ~80 ks between the observations 
(ObsID 7290 and 5436). 

We reprocessed these datasets with CIAO 4.1 and com- 
bined the 2 datasets with the script merge-all. We then 
used wavdetect to detect sources in the energy ranges 

4 http://sirnbad.u-strasbg.fr/ 
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Figure 7. Count rate vs. hardness ratio for the Chandra sources 
detected in the M13 field of view. For clarity, only the 30 X- 
ray sources closest to the center are labelled. Hardness ratios are 
computed from the counts in the bands Bl = 0.3 — 1.5 keV and 
B2 = 1.5 — 8 keV. Black solid lines indicate the hardness ratios 
corresponding to a power law model (PO) with different photon 
indices and a bremsstrahlung emission (BR) with different tem- 
peratures kT in keV. The luminosity for a source in the cluster 
with a spectrum modeled by a power law of photon index 1.5 is 
indicated on the right axis. The mean errors for 3 luminosity levels 
are indicated by gray bars. 

0.5 — 2 and a broader 0.3 — 8 keV band with scales of 1.0, 
1.4, 2.0, 2.8, 4.0, and 5.6 pixels. We selected a thresh- 
old probability of 10 -6 , designed to give one false source 
per 10 6 pixels. We found 47 X-ray sources in a 9' x 9' 
region centered on M13 (Table H]). Two sources are only 
detected in the energy range 0.5 — 2 keV (X14 and X30) 
and 12 are only detected in the broad band (including 
X3, X10 and X13 inside the half- mass radius). This list 
includes the 14 sources dete cted previously with XMM- 
Newton (Gendre e t al.l 



catalog (Wa tson et al. 



2003) and available in the 2XMM 



20091 ). An image that covers the 



half-mass radius of M13 is shown in Figure [6l A hardness 
ratio diagram is given in Figure [71 

In the energy range 0.5 — 2 keV, a total of 15 sources 
are located inside the half-mass radius, and 8 of them are 
inside the core radius. Within a circle with a radius of 
4.'2, we detected 35 sour c es (0.5 — 2 keV band). Using the 
study of Hasinger et al.l ([2005), which gives an estimate 
of the space density of background sources (active galac- 
tic nuclei, AGN, in the soft X-ray band 0.5 — 2 keV), 
we estimated the expected number of AGN for several 
regions (see Table [3|). We assumed that the minimum 
flux of a detected source during the observation corre- 
sponds to 5 counts. A clear excess of detections is ob- 
served in the direction of M13. We limited this study to 
a region of radius 4 .'2 where the sensitivity drops to 86% 
of the central sensitivity and the PSF is slightly more 
elongated. As we assumed a constant minimum flux in 
this region, this could have led to a possible deficit of 
detections compared to the expected number of back- 
ground sources in the larger region. However, this deficit 
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Table 4 

List of Chandra X-ray sources in the field of M13 





RA 


Dec 


Q^vr> prrnr 
u <j /(j ci i ui 


J_-/±k5 lid-Llv^C 


Npf Clnnnt Rafp 

IMCli V^U Lllllj ±V,CXLC 


R1 C!nnnt<3 

J_J J_ VUU11UB 


R9 C!nnnt<3 


HR 




[h;m:s] 


[d:m:s] 


P'l 




ryin- 3 r.-ni- q — ll 

[ x iu circ s J 


u . o i . o Ke v 


i . o o Ke v 




XI 


16:41:42.015 


36:27:34.65 


0.70 


0.13 


0.24±0.07 


10 


5 


-0.33±0.36 


X2 


16:41:42.676 


36:27:47.65 


0.68 


0.24 


0.23±0.07 


11 


4 


-0.47±0.35 


X3 


16:41:40.366 


36:27:42.84 


0.68 


0.26 


0.12±0.05 


4 


3 


-0.14±0.53 


X4 


16:41:42.401 


36:27:28.29 


0.50 


0.26 


0.10±0.04 


6 





-1.00±0.41 


X5 


16:41:43.160 


36:27:25.44 


0.42 


0.40 


0.98±0.14 


22 


37 


0.25±0.18 


X6 


16:41:42.480 


36:28:07.52 


0.15 


0.48 


3.19±0.24 


99 


77 


-0.12±0.11 


X7 


16:41:43.761 


36:27:57.84 


0.08 


0.51 


11.05±0.45 


575 


33 


-0.89±0.05 


X8 


16:41:41.078 


36:27:08.11 


0.69 


0.56 


0.14±0.05 


6 


6 


0.00±0.41 


X9 


16:41:44.822 


36:27:37.04 


0.26 


0.64 


1.66±0.18 


46 


47 


0.01±0.15 


X10 


16:41:38.839 


36:28:04.53 


0.40 


0.69 


0.05±0.03 


1 


2 


0.33±0.80 


Xll 


16:41:44.622 


36:28:15.97 


0.39 


0.84 


0.23±0.07 


12 


1 


-0.85±0.34 


X12 


16:41:38.439 


36:28:23.48 


0.55 


0.96 


0.36±0.08 


14 


8 


-0.27±0.30 


X13 


16:41:46.283 


36:28:01.74 


0.69 


1.00 


0.06±0.04 


2 


4 


0.33±0.57 


X14 


16:41:42.775 


36:28:40.55 


1.36 


1.02 


0.13±0.05 


2 


5 


0.43±0.52 


X15 


16:41:46.838 


36:27:30.93 


0.31 


1.06 


0.77±0.12 


14 


30 


0.36±0.21 


X16 


16:41:39.393 


36:26:35.67 


0.58 


1.17 


0.27±0.07 


9 


7 


-0.12±0.35 


X17 


16:41:39.950 


36:26:28.15 


0.31 


1.26 


1.19±0.15 


51 


17 


-0.50±0.17 


X18 


16:41:38.209 


36:26:26.71 


0.12 


1.41 


6.39±0.34 


228 


124 


-0.30±0.07 


X19 


16:41:33.832 


36:28:14.33 


0.41 


1.67 


0.57±0.10 


19 


13 


-0.19±0.25 


X20 


16:41:34.074 


36:26:58.57 


0.84 


1.67 


0.13±0.05 


6 


2 


-0.50±0.48 


X21 


16:41:49.019 


36:26:44.73 


0.65 


1.75 


0.27±0.07 


7 


10 


0.18±0.34 


X22 


16:41:42.346 


36:25:49.57 


0.09 


1.86 


10.76±0.45 


460 


124 


-0.58±0.06 


X23 


16:41:45.407 


36:29:35.47 


0.52 


2.06 


0.05±0.03 


2 


1 


-0.33±0.80 


X24 


16:41:49.794 


36:26:24.47 


0.62 


2.08 


0.20±0.06 


7 


5 


-0.17±0.41 


X25 


16:41:41.926 


36:25:26.92 


0.62 


2.23 


0.50±0.10 


12 


19 


0.23±0.25 


X26 


16:41:45.947 


36:25:32.48 


0.68 


2.31 


0.27±0.07 


14 


5 


-0.47±0.31 


X27 


16:41:50.006 


36:25:42.67 


0.65 


2.59 


0.14±0.05 


6 


2 


-0.50±0.48 


X28 


16:41:37.792 


36:25:05.22 


0.20 


2.70 


2.84±0.23 


90 


69 


-0.13±0.11 


X29 


16:41:27.839 


36:27:34.26 


0.90 


2.78 


0.06±0.04 


1 


3 


0.50±0.68 


X30 


16:41:41.646 


36:24:41.70 


1.44 


2.98 


0.08±0.04 


4 


1 


-0.60±0.60 


X31 


16:41:50.285 


36:25:13.85 


0.25 


3.00 


1.30±0.16 


62 


8 


-0.77±0.15 


X32 


16:41:54.727 


36:26:04.15 


0.75 


3.09 


0.43±0.09 


8 


18 


0.38±0.27 


X33 


16:41:45.948 


36:24:22.01 


0.62 


3.42 


0.69±0.12 


33 


14 


-0.40±0.20 


X34 


16:41:48.727 


36:30:48.17 


0.75 


3.43 


0.41±0.09 


13 


13 


0.00±0.28 


X35 


16:41:32.606 


36:30:37.58 


0.37 


3.46 


2.17±0.20 


33 


93 


0.48±0.12 


X36 


16:41:36.849 


36:31:15.09 


1.08 


3.70 


0.28±0.08 


9 


10 


0.05±0.32 


X37 


16:41:39.239 


36:23:43.84 


0.62 


3.98 


1.14±0.15 


51 


21 


-0.42±0.16 


X38 


16:41:30.540 


36:31:02.52 


1.59 


4.03 


0.16±0.06 





10 


1.00±0.32 


X39 


16:41:39.846 


36:23:31.32 


1.20 


4.17 


0.21±0.07 


12 


7 


-0.26±0.32 


X40 


16:41:21.629 


36:29:01.59 


1.35 


4.24 


0.10±0.05 


1 


7 


0.75±0.46 


X41 


16:42:02.732 


36:27:21.40 


1.09 


4.25 


0.14±0.05 


5 


8 


0.23±0.39 


X42 


16:41:52.258 


36:23:23.67 


1.13 


4.79 


0.11±0.05 


3 


5 


0.25±0.50 


X43 


16:41:57.158 


36:31:24.07 


1.57 


4.86 


0.26±0.08 


6 


11 


0.29±0.34 


X44 


16:41:32.745 


36:23:08.08 


1.86 


4.88 


0.16±0.06 


8 


21 


0.45±0.26 


X45 


16:42:03.911 


36:25:27.95 


1.60 


5.00 


0.10±0.04 


7 


1 


-0.75±0.46 


X46 


16:42:02.659 


36:30:59.96 


1.17 


5.37 


0.35±0.08 


11 


13 


0.08±0.29 


X47 


16:42:03.380 


36:24:15.48 


0.42 


5.55 


4.93±0.31 


126 


174 


0.16±0.08 



Note. — For each source we give its ID, coordinates and position error (columns 1 to 4). We report the distance to the center in 
column 5. We then give the net count rate, the number of counts in the 2 bands Bl = 0.3 — 1.5 keV and B2 = 1.5 — 8 keV, and the 
hardness ratio HR = (B2 Bl) / (Bl + B2). 

120091 ) gives a flux in the band 0.2-12.0 keV of 3.9 ± 0.7 x 
10 -14 erg s _1 cm -2 (8.8 x 10 31 erg s _1 at the distance 
ofM13). 

The position of the corresponding Chandra source (X6, 
see gSJ) is RA = 16Ml m 42 s .48, Dec = +36°28 / 07 / /6 
(J2000). We comput e d the 95% error from the formula 
given by iHong et all (j2005f ) and found / /32. Including 
the absolute pointing error of Chandra ACIS-S (Off 7 at 
95%) we find a total error on the position of f !77 at 95% 
confidence level. The variable star 4 lies 1//15 away with 
a positional error of l'/2 at 95%. 

We extracted the lightcurve of source X6 from the 
Chandra data using the CIAO tool dmextract (see Fig- 
ure [8j). The count rate of source X6 was higher by a 
factor ~3 at the beginning of the observation and for 
~10 ks. The background is negligible throughout the 
observation. 



is not observed, strengthening the excess of detections 
observed which may thus be slightly underestimated. As 
Ml 3 is at a Galactic Latitude of 40.91° we expect only 
few foreground X-ray sources. Therefore, most of the X- 
ray sources inside the half-mass radius are likely to be as- 
sociated with M13, and a tentative estimate of 19.6 ±3.9 
sources could be associated with the cluster. 

4. DETAILED RESULTS FOR THE VARIABLE STAR 4 

4.1. X-ray counterpart: lightcurve and spectrum 

The star 4 falls in the error radius of an X-ray source 
det ected with XMM-N ewton on 2002 January 28 and 
30 (jGendre et al.ll2003h . The X -ray source is described 
as varia ble by at least a f actor ~2 on timescale of 
years by iGendre et al.l (j2QQ3h as it should have been de- 
tected ^_by_D^evious X-ray observations wi th ROSAT in 
1994 (IVerbuntll200Th . The 2XMM catalog (IWatson et al.l 
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v37 




Figure 8. X-ray lightcurve of Chandra source X6 in M13. The 
lightcurves start time is 2006-03-09 23:01:13 UT (MJD 53803). The 
mean level at low luminosity is indicated by a dashed line. 



Spectrum of Chandra source X6 in M13 




Figure 9. X-ray spectrum of Chandra source X6 in M13. The 
solid lines are the absorbed bremsstrahlung model. The residuals 
are indicated in the bottom panel 

We extracted the spectra of X6 for both observations 
using the tool psextract. We generated the RMF re- 
sponse files separately with the tool mkacisrmf. We 
used Sherpa to fit both spectra simultaneously, exclud- 
ing the first 10 ks of data. Given the limited number 
of counts (a total of 134 counts), we tried simple mod- 
els with a fixed absorption of 1.1 x 10 20 cm -2 to fit the 
spect rum. We used the TBabs model with abundances 
from I Wilms etall (pQQOh . Good fits were obtained for a 
power law with an index of 1.3±0.3 (x 2 /dof = 3.9/7) or a 
bremsstrahlung model with temperature of kT ~35 keV 
(x 2 /dof = 4.2/7). The temperature is however not well 
constrained (kT < 50 keV at la). The spectrum and 
best fit are reported in Figure [9] 

Based on the power-law models, we estimated 
that the flux of the source in the band 0.3- 
8.0 keV (with la error) for the second observation 
is 2.2 ±0.8 x 10" 14 erg s _1 cm -2 (5.0 x 10 31 erg s _1 
at the distance of M13). During the increase of 
counts in the first observations, the flux reached 



cm 2 (1.1 x 10 32 erg s 



at 



5.2 ± 1.9 x 10- 14 erg s 
the distance of M13). 

Due to the low number of counts, we estimated the 
median energy at 3 different tim e intervals to lo ok for 
possible spectral variability (see iHong et al]|2004 for er- 
ror estimations). During the first ~ 10 ks, the median 
energy is 1.29 ±0.12 keV, then it rises to 1.77 ±0.25 keV 
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Figure 10. Close up Faulkes U-band image around the Chandra 
source X6 and the variable star 4. Variable star candidates and 
X-ray sources are indicated by their 95% error radius (black and 
red respectively), and known variables are indicated with a blue 
dashed circle of radius 1//8 for clarity. Left: reference image. Right: 
combination of subtracted images as presented in the text. We also 
discuss in the text complex features seen in the background. 

at the end of the first observation, possibly indicating 
a slight hardening of the source in the low luminosity 
state. In the second observation the median energy is 
1.26 ± 0.16 keV. 

4.2. Maximum U-band magnitude 

We averaged the subtracted images taken between 
MJD 53419 and 53433 with the FTN, when the star 
was at a maximum of luminosity. We also averaged the 
images between MJD 53472 and 53477 during the mini- 
mum. In those two images, the source appears as positive 
and negative, respectively, compared to the emission of 
the source in the reference image. In order to obtain 
the total flux of star 4 between the minimum and maxi- 
mum states, we subtracted the negative image from the 
positive image. In Figure [lOj we show the reference im- 
age and the obtained combination of subtracted images 
around star 4. An excess is clearly seen, and the error 
circle of the variable star and the X-ray source overlap, 
indicating that the two sources are counterparts. One 
can note that the background is irregular, with many 
features that show positive and negative peaks. Those 
are residuals from the subtraction of stars. They can 
be ruled out as variable stars as the integrated flux in 
a sufficiently large circle generally remains close to zero 
and well below the standard deviation of the surround- 
ing background. The shape of the emission is also not as 
point-like as for star 4, and does not have the expected 
PSF size. 

We studied the stack of subtracted images in the U- 
band with IRAF imexamine in order to determine the 
magnitude of the star. We first set the zero point with 
the reference image. We used 6 isolated and bright ref- 
erence stars with a known spectral type from SIMBAD. 
These are not classified as variable and possess B and 
V magnitude measurements. Based on this informa- 
tion, we extrapolated the expected U magnitude of those 
stars, and scaled the measured U magnitude accordingly. 
The resulting standard deviation for the set of reference 
stars is <ttj = 0.3. As the extinction is low in this field, 
E(B — V) = 0.02, it is unlikely to have an effect on this 
measurement. 

Based on the scatter of the background in one sub- 
tracted image, the size of the FWHM and the zero point, 
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Table 5 

Hubble Space Telescope archived observations. 



Prop. ID Date MJD Filter Exp. time 



10349 2005-08-15 53597 F435W (B) 


2 x 680s, 


1x120s 


F625W (R) 


1x360s, 


1x20s 


F658N (Ha) 


1x800s, 


1x690s 


10775 2006-04-02 53827 F606W (R) 


1x567s 




F814W (I) 


1x567s 





we estimate that the 5<r detection threshold corresponds 
to U = 19 ± 0.3. We derived a U magnitude of 17.3 ± 0.3 
for star 4. The star has not been reported before, so it is 
likely that the star was fainter and below the detection 
threshold during its minimum of luminosity. The mini- 
mum flux would then be more than 1.7 magnitudes lower 
than the maximum flux. Given the large error bars, we 
can assume that the measured magnitude in the combi- 
nation of subtracted images can be directly taken as the 
maximum magnitude reached by the star. The absolute 
magnitude is Mtj = 3.0 ± 0.3 if the object is associated 
with M13. 

4.3. Search for a counterpart with Hubble 

M13 has been observed with the Hubble Space Tele- 
scope for 6 programs after 1995. We checked the previews 
of the images acquired and found 3 programs (proposal 
IDs 10775, 10349, and 8174) where the region of the X- 
ray source X6 is covered. The dataset with proposal ID 
8174 is composed of short exposures only (e.g. 26 s) so 
we focused on the 2 other datasets obtained with the 
Advanced Camera for Surveys (ACS) and its Wide Field 
Camera (WFC) detector. The list of available exposures 
is given in Table [5] 

We downloaded the pipeline-reduced images (drz files) 
and corrected the astrometry with GAIA using 3 refer- 
ence stars in the field of view and out of the crowded 
region of the cluster from the UCAC2 catalog, giving an 
RMS of 0'/0003. The error thus mainly comes from the 
UCAC2 reference stars, which is 0^04 (at la) in this field 
of view. Images of the region around star 4 are given in 
Figure [H 

We detected the stars and performed the photometry 
measurement with IRAF/DAOPHOT. We used 6 iso- 
lated, non-saturated stars to create the PSF model for 
each image, and applied the task allstar to fit the stars 
per group of at most 60 stars. The photometry was cal- 
ibrated using the zero points provided on the ACS web- 
pag4E to get Vega magnitudes in each filter. We detected 
11 counterparts inside the Chandra error circle of X6 in 
the F814W image, and only 5 in the F658N image. We 
obtained magnitudes in the 4 broad bands for 7 stars, la- 
belled SI to S7 according to their distance to X6, which 
are indicated in Figure [T21 They are also reported in the 
3 different color-magnitude diagrams (CMDs) shown in 
Figure [TTJ Each diagram was generated with 2 images 
taken during the same program. 

We note that some spikes from bright stars are ap- 
parent in the images (see Figure [T2]) . Unfortunately, 
these spikes likely contaminate SI, S2 and possibly S4 
in the first dataset (F625W and F658N), and S4, possi- 
bly SI in the second (F606W and F814W). This may 

5 http:/ /www. stsci.edu/hst/acs/analysis/zeropoints 



explain the extreme red color of SI and S2 in the 
CMD (F625W, F435W-F625W) which are then bluer 
and fainter relative to the main sequence in the CMD 
(F606W, F606W-F814W). Maybe the blue color of S4 in 
the CMD (F625W, F435W-F625W) is also due to con- 
tamination as the star is faint. There remains, however, 
a possibility that SI, S2 and S4 showed color variability. 

The CMD (F625W, F658N-F625W) can be used as an 
indicator of Ha excess. We note that compared to the 
distribution, SI seems to present a ~3cr Ka excess. There 
does not seem to be a strong contamination in the F658N 
band, and it is likely that it would be compensated or 
dominated by the contamination in the F625W band. 

5. DISCUSSION 

5.1. Star 4- & DN outburst 

Using observations from different observatories in X- 
ray and optical, we find compelling evidence that the 
variable star 4 (X-ray source X6) is a DN. First, the 
star shows variability on a timescale of months, con- 
sistent with DN outbursts that could last at least 24 
days and reached a maximum absolute U magnitude of 
Mtj = 3.0 ± 0.3 if it is located in M13. We do not 
have color information, but we note that known DNe 
i n outburst can hav e U-V values between and -1 
(jHarrison et a l. 2004). Therefore star 4 is likely to have 
a maximum magnitude of My > 3.0. Previously studied 
DNe have a maximum absolute ma gnitu de in outburst 
of M v ~ 2.7 - 5.2 (jHarrison et al.M2QQ4f ). The optical 
brightenings of the variable star 4 would thus be consis- 
tent with bright outbursts from a DN located in the core 
of M13. 

The position of this star is consistent with an X-ray 
source detected using ROSAT, XMM-Newton and Chan- 
dra. The X-ray source has varied in flux between the 
ROSAT and XMM-Newton observations as expected for 
CVs, but many other X-ray sources show variability. The 
emission of the Chandra X-ray source is compatible with 
the expected emission of a DN. Indeed, DN have lumi- 
nosities in the range 10 29-33 erg s -1 and bremsstrahlung 
temperatures of kT - 10 keV (e.g. IVerbunt et all 119971 : 
Bask ill et al.ll2005l ). It is thus similar in luminosity to 
the brightest DN such as U Gem or SU UMa (e.g. 
Bask ill et al.ll2005l ). We note that the temperature we 
estimated is also consistent with the higher tempera- 
tures found^Jor_inter mediate polars, e.g. 20-30 keV 
(e.g. Ezuka fc Ishidal fl999h . A variation in the X-ray 
count rate of a factor 3 has been observed with Chan- 
dra. This kind of fluctuation is common for CVs or 
magnetically active binaries. However the latter gener- 
ally show an exponential decay of the X-ray flux over 
a time range of typically 10 ks which is not seen here. 
We note that this X-ray drop occurred 2 days before U- 
band observations with the FTN which show the star 
in a bright state. This behavior seems consistent with 
the X-ray /optical anticorrelation observed for some DNe. 
For example, SU UMa showed a sudden drop in X-rays 
by a factor 4 about half a day after the optical rise 
flCollins fc Wh eatlev 2010). The hard X -ray emission of 
SS Cvg (|McGowan et al.l l2004) is almost suppressed dur- 
ing the optical outburst. This is somewhat contrary to 
the possible hardening of the source that we observed 
after the drop in flux. However, SS Cyg also showed an 
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Figure 11. Color magnitude diagrams of M13 from HST ACS/WFC images. When detected, the stars SI to S7 shown in Figure fT2l are 
indicated. Errors bars are shown on the left of the diagram, which correspond to the mean error for stars in the corresponding magnitude 
range. In the first diagram, the mean, la and 3a envelopes of the distribution are shown as dashed lines. 



intermediate state where si milar hardening and dec rease 
of the X-ray flux occurred (|McGowan et al.l 12004). SS 
Cyg also showed an X-ray flux increase at the beginning 
of the optical outburst for a bout half a day, as exp ected 
from boundary-layer models (Wh eatlev et al.ll2003[ ). It is 
possible that the drop in flux we observed corresponds in 
fact to this increase that precedes the optical outburst. 

Finally, using archived data from the Hubble Space 
Telescope we found ~11 stars inside the Chandra er- 
ror circle. No bright star was found, thus the object 
is likely to be in quiescence during those observations. 
One of the brightest counterparts (SI) presents a possi- 
ble Ha excess. This kind of emission could come from 




2005-08-15 / MJD 53597 2006-04-02 / MJD 53827 



Figure 12. HST ACS/WFC snapshots. Images from the datasets 
with proposal ID 10349 and 10775 are on the left and right columns 
respectively. The size of each image is 4'/ 2 x 2 / /5, North is up, East 
left. The 95% Chandra error on the position of the candidate CV 
are indicated as red circles (small dashed: only detection error, 
large plain: including absolute pointing error). SI to S7 are shown 
with a 2 pixels radius circle. 



an accretion disk around a white dwarf and would indi- 
cate that SI is the true counterpart. However, SI does 
not show a blue excess as it is generally expected from 
DNe, but we note that some DNe in quiescence in GCs 
are c onsistent with main sequ ence stars (e.g. CV2 in 
M22. iPietrukowicz et all 120051 ). We note that the star 
S4 is suspiciously blue in one observation and could have 
shown color variability, however, it was close to a spike 
from a bright star and measurements may not be reli- 
able. We could measure magnitudes and colors for the 
7 brightest counterparts only. Thus, one of the fainter 
counterparts, for which we could not determine reliable 
colors, could turn out to be a more convincing associa- 
tion. 

If we consider the most studied DNe, their absolute 
magnitude in qu iescence lie in the range My = 6—10 (e.g. 
lPattersonll2011[ ) . The amplitude of the outburst is gener- 
ally dM v =2 — 6. DNe with brighter outbursts also tend 
to have lower amplitude. At the distance of M13, the qui- 
escent magnitude range would be V = 20.2 — 24.2. This 
outburst being bright, we would then expect a counter- 
part in quiescence to have a V magnitude in the brighter 
part of this range. This makes it less likely for fainter 
stars such as S4 to be the true counterpart, while SI 
stands as a good candidate. However, further observa- 
tions are still needed to locate the true counterpart, with 
better signal to noise ratio and less contamination from 
bright stars in order to confirm an Ha excess. 

5.2. Population of CVs in Ml 3 

Based on the encounter rate in M13 and the correla- 
tion of this parameter with the nu mber of X-ray sources 
reported bv iPoolev fc Hutl ([2006) , we expect about one 
source within the half- mass radius of Ml 3 above a lumi- 
nosity threshold of 4 x 10 31 erg s _1 , where cluster X-ray 
sources are dominated by CVs. Inside the half-mass ra- 
dius, we found 2 X-ray sources above this threshold. One 
of them is very so ft and already repor ted as a quiescent 
LMXB candidate (jGendre et all 12003( 1. The other is the 
DN discovered in this work. This is consistent with the 
estimation and strengthens the identification of X6 as a 
CV. 
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llvanova et al.l J2006) performed simulations based on 
different CV formation channels. Their low-density clus- 
ter is the closest to M13, and they estimated that a total 
of 156 CVs can be formed in such a cluster. From this 
simulation, about 40% of the CVs formed are purely pri- 
mordial CVs (i.e. no dynamical interaction were needed 
to lead to their formation). Another 46% of the CVs are 
also primordial, but entered the core of the cluster be- 
fore the common envelope phase. Only 5% of the CVs 
in the simulation come from formation channels that in- 
volve collisions or binary exchange mechanisms. There- 
fore, the population of CVs in M13 could more closely 
resemble the population of field CVs than populations 
in denser GCs, where dynamical effects are believed to 
modify and enhance the production of CVs. 

There have been few detections of DNe in GCs, with 
only 13 DNe reported so f a r in the 157 Galactic GCs (e.g . 
iPietrukowicz et al.l 120081 : iKaluznv fc Thompson! l2QQ9h . 
Based on our knowledge of field CVs, which is extremely 
biased, it is tempting to try to estimate the expected 
number of DN outbur sts we should have d etected in Ml 3, 
for example following Sh ara et all (jl996[ ). The DN pop- 
ulation constitute s about half of the k nown CVs listed 
in the catalog of iDownes et al.l (j2QQlh . which contain 
mostly field CVs. We could then expect ~78 DNe in 
M13, based on those rough estimations. The mean duty 
cycle, or probability of a DN being in outburst at a ran- 
dom epoch, has been found to be ^15% base d on the 
observation of 21 DNe (jSzkodv fe Matteil fT984). We can 
group our observations into 5 epochs of few days, thus 
the likelihood of a DN being in quiescence in those 5 
epochs is 0.85 5 = 0.44. Thus -42 DNe should have 
been in outburst in at least one of the epochs. It is 
then harder to estimate the number of outbursts that 
would have been bright enough to be detected in our ob- 
servation. Our detection limit at U ~ 19 would convert 
to V ~ 20 (M y — 5.7) due to the g eneral blue color of 
DN outbursts (Harri son et al.ll2004f ). We can then use 
the distribution of V magnitudes at maximum light from 
iPattersonl (|2Qllh based on the observations of —45 DNe, 
which shows that we could have expected to detect at 
least two third of those DNe in outburst at the distance 
of M13. Finally, — 28 DNe in outburst would have been 
expected when we detect only one. 

However, it is important to note that this estimate is 
based on several approximations, which makes this num- 
ber likely to be well over-estimated. If we suppose that 
the mean duty cycle of DNe is 1% instead of 15%, we 
then expect —2 DNe in outburst. We can also suppose 
that the proportion of DNe is much less than half of all 
CVs. It would have to be less than 5% to expect so few 
DNe in eruption, as in our observations. Considering 
that M13 is not an extremely dense GC, with most of 
its CVs likely being primordial, we may be able to de- 
duce basic information on the composition of the global 
Galactic CV population. Even if there are large errors 
in each of our estimates, it seems likely that the mean 
duty cycle of DNe is lower than 15%, and/or the intrinsic 
fraction of DNe among CVs may be smaller tha n cata - 
logs suggest. If the predictions of llvanova et al.l {2006) 
regarding total numbers of CVs in GCs are correct, then 
a possible solution to the apparent discrepancy between 
M13 and the field could be that there is a sub-population 
of CVs that undergo outbursts but have lower luminosi- 



ties and longer outburst recurrence times than those of 
the currently known field DNe. 

Several population synthesis models already showed 
that a sub-population of faint CVs with long outburst 
recur rence times may dominate the total number o f 
CVs (|Kolblll993l : lHowell et al .1119971 : 1 Willems et al.ll20Q5l ). 
However there are still few e xamples of such a sub- 
population (e.g. Gans icke et al.ll2009l ). Deeper surveys of 
the Galactic Plane are thus needed to catch fainter CVs 
in order to test this hypothesis. Having observations on 
timescales of decades could also help to detect outbursts 
with longer recurrence times. This will maybe be possi- 
ble through the DASCH project which aims to digitize 
half a million photogr aphic plates obtained over 100 years 
in the last century (iGrindlav et al.l 120091 : lLavcock et al.l 
l2010l : IServillat et al.ll201lI T^ 

On the other hand, it is possible that our observations 
were inefficient in detecting DN outbursts. Considering 
the expected recurrence time of months to years and a 
duration of at maximum a few weeks for DN outbursts, 
the optical observation coverage that is likely to detect 
most of the DN outbursts requires both short scale ob- 
servations (over few nights) and a recurrence every few 
weeks, if possible over a year. This kind of survey is thus 
time consuming, particularly for the large telescopes that 
are needed if the target is crowded as it is the case for 
GCs. 

Future variability surveys, such as the VISTA Vari- 
ables in The Via Lactea (VVV, iMinniti et all I2010D 
public s urvey and the La rge Synoptic Survey Telescope 
(LSST, E vezic et al. 2008) will conduct an even denser 
coverage with a 4m and 8.4m telescopes, respectively. 
One objective would be to uncover more DNe in GCs 
and in the field, in order to base our study of CV popu- 
lations on samples with a decent number of objects. 

6. SUMMARY AND CONCLUSIONS 

We performed a search for optical outbursts in the GC 
M13. We found one star in M13 which showed vari- 
ability consistent with outbursts. This star is consistent 
with being a DN, based on the optical variability and 
flux, and on the X-ray variability, spectrum and lumi- 
nosity. Several possible counterparts were detected in 
Hubble images, including one with a possible Ha excess 
and consistent magnitude in quiescence. 

The fact that we detect one DN when more could have 
been expected suggests that the known population of lo- 
cal CVs is biased, and short-period, faint, low- duty- cycle 
(1%) DNe are likely to dominate the global CV pop- 
ulation, as expected from population synthesis models. 
Their proportion is however not well constrained and fu- 
ture variability surveys will help to better characterize 
the global CV population. 
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